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SUMMARY 


The  elastic  constants  of  a  fiberglass  epoxy  unidirectional 
composite  are  determined  by  measuring  the  phase  velocities  of 
longitudinal  and  shear  stress  waves  via  the  through  transmission 
ultrasonic  technique  .  The  waves  introduced  into  the  composite 
specimens  were  generated  by ypiezoceramic  transducers.  Geometric 
lengths  and  the  times  required  to  travel  those  lengths  were  used 
to  calculate  the  phase  velocities.  The  model  of  the  transversely 
isotropic  medium  was  adopted  to  relate  the  velocities  and 
elastic  constants. 
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INTRODUCTION 


Filamentary  composites  have  proved  their  usefulness  in 
modern  aircraft,  spacecraft,  land  and  water  vehicles,  and  many 
other  structural  applications. 

One  of  the  special  characteristics  of  these  materials  is 
their  anisotropic  elastic  behavior,  that  is,  their  elastic 
properties  which  change  with  orientation.  Knowledge  of  such 
properties  is  essential  both  for  optimal  design  and  for 
evaluation  of  performance.  Elastic  properties  are  usually 
represented  by  a  set  of  constants  which  can  be  determined  by 
conventional  stress-strain  measurements  in  different  loading 
configurations  (as,  for  instance,  in  standard  tensile  tests)  or 
by  ultrasonic  methods. 

Ultrasonic  testing  is  a  nondestructive  evaluation  (NDE) 
technique  in  which  the  propagation  characteristics  of  stress 
waves  in  a  material  are  determined  .  The  wave  transmission 
characteristics  may  then  be  related  to  such  factors  as  the 
elastic  properties  of  the  material  or  the  presence  of  flaws. 
Ultrasonic  methods  have  been  extensively  used  for  the 
determination  of  velocities  of  stress  waves,  elastic  constants  of 
the  stiffness  matrix,  fiber  fraction,  degree  of  fiber 
orientation,  presence  of  voids,  and  attenuation  and  damping 
parameters  [1-6].  One  of  the  main  advantages  of  ultrasonic 

t 

methods,  as  they  are  essentially  nondestructive,  is  their 
suitability  for  the  evaluation  of  ”in  service”  components  . 

The  purpose  of  this  work  is  to  establish  the  experimental 
values  of  the  velocities  of  ultrasonic  waves  of  a  specific 


1 


composite  laminate  of  fiberglass  epoxy  in  which  the  fibers  are 


continuous 
the  values 
determined . 
purposes 
performance 


and  unidirectional.  From  the  measured  velocities, 
of  elastic  constants  of  the  stiffness  matrix  are 
Such  parameters  can  then  be  used  both  for  design 
or  as  standards  for  the  NDE  evaluation  of  the 
of  components. 
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TRANSVERSELY  ISOTROPIC  COMPOSITE 

In  filamentary  composite  materials,  the  presence  of  a 
distinctive  fiber  orientation  gives  rise  to  elastic  anisotropy. 
Unidirectional  fiber  composites  consist  of  a  parallel  set  of 
cylindrical  fibers  embedded  in  an  otherwise  homogeneous  matrix 
material.  If  the  constituent  materials  are  isotropic,  the 
resulting  composite  will  be  classified  as  transversely  isotropic, 
or  hexagonally  symmetric. 

A  transversely  isotropic  material  can  be  formed  from 
prepreg  sheets  of  continuous  parallel  fibers  that  are  layered  to 
form  a  unidirectional  composite  [7] .  A  schematic  of  this 
configuration  is  shown  in  Fig.1.  When  the  transversely  isotropic 
model  is  adopted,  the  fibers  are  assumed  to  be  randomly  arranged 
in  the  xy  plane. 

To  characterize  the  mechanical  behavior  of  transversely 
isotropic  materials,  five  independent  elastic  constants  are 
required  [8]  . 
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BASIC  PRINCIPLES 


The  elastic  properties  of  a  filamentary  transversely 
isotropic  composite  are  evaluated  by  measuring  the  phase 
velocities  of  longitudinal  and  shear  plane  waves  propagating 
through  the  material  at  various  angles  with  respect  to  the 
fiber  direction. 

Depending  on  the  direction  cosines  chosen  for  the  normal  of 
the  propagating  wave,  the  equations  yielding  the  five  elastic 
constants  in  terms  of  velocities  can  be  obtained.  These  equations 
have  been  derived  for  transversely  isotropic  materials  [9] •  For 
a  cartesian  coordinate  system  x,  y  and  z  in  which  2  coincides 
with  the  fiber  direction  in  the  transversely  isotropic  material, 
the  relations  between  the  elastic  constants  of  the  stiffness 
matrix  and  phase  velocities  are  summarized  in  Table  1 .  The 
propagation  of  waves  in  any  direction  or  polarization  direction 
of  a  wave  is  represented  by  its  direction  cosines  with  respect  to 
the  X,  y,  and  z  directions  (for  example  (1,0,0)  defines  the  x 
direction).  P  denotes  an  imposed  longitudinal  wave,  for  which 
propagation  and  polarization  directions  are  coincident;  SH 
denotes  an  imposed  shear  wave  with  polarization  in  the  y  axis  ; 
and  SV  denotes  imposed  shear  wave  with  polarizaton  in  the  x-z 
plane . 

The  constants  C,,,  and  can  be  determined  by 

11  33  44  00 

direct  measurement  of  the  velocities  of  appropriate  waves 
propagating  in  the  principal  directions  x,y  and  z.  The 
remaining  constant,  »  can  be  determined  for  a  wave 
propagating  along  an  intermediate  direction  between  the  axes  x 
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and  z  .  (Table  1  shows  the  equations  for  some  intermediate 

directions,  defined  by  their  angles  with  respect  to  the  z  axis.) 

The  constant  C.-  can  be  calculated  from  C..  and  as 

12  11  66 

as  2Cgg  . 

The  model  adopted  describes  the  composite  as  homogeneous,  so 
the  density  is  assumed  to  be  constant.  The  model  is  valid  for 
stress  waves  having  wavelengths  that  are  large  compared  with  the 
thickness  of  the  prepreg  sheets.  Previous  studies  [10,11 ]  have 
shown  that  waves  can  be  scattered  by  reflection  if  the 

wavelengths  have  dimensions  comparable  to  the  thickness  of  the 

lamina  or  the  fiber  diameter.  Effects  of  dispersion  may  be 
present  in  such  cases,  and  a  better  model  ought  be  a  layered 
medium  model  instead  of  the  homogeneous  model. 

For  the  homogeneous  model,  the  correspondence  between  the 
elastic  constants  of  the  stiffness  matrix  and  the  magnitudes  of 
extenslonal  moduli,  shear  moduli,  and  Poisson's  ratios  are  shown 
in  Table  2. 

For  composites  with  symmetries  other  than  hexagonal,  the 
procedure  of  measuring  velocities  can  also  be  used,  but  in  these 
cases  a  different  set  and  number  of  independent  elastic  constants 
must  be  determined  .  It  would  be  possible  to  determine  the 
symmetry  axes  (if  any  )  of  any  specimen  by  making  sufficient 
measurements.  This  would  involve  a  large  amount  of  experimental 
work  and  analysis  which  could  be  greatly  reduced  if  some  prior 
information  is  known  about  the  symmetries  of  the  material. 
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EXPERIMENTAL  PROCEDURES 


Procedures  and  Equipment 

The  measurements  were  made  by  the  through  transmission 
technique  using  specimens  of  different  lengths  for  each  of  the 
measurement  directions.  A  schematic  of  the  measuring  system  is 
shown  in  Fig. 2  .  The  system  consisted  of  a  pulsed  oscillator 
(Arenberg  model  PG-652)  for  generating  sinusoidal  waves  ,  a  low 
frequency  inductor  (Aremberg  model  LFT-500  ),two  broadband  (0.1 
to  3*0  MHz)  transducers  (AET  model  FC-500)  for  longitudinal  waves 
with  approximately  flat  sensitivity  of  -85  dB  (relative  to  1V/ 
Bar),  two  transducers  (Panametrics  model  V154)  shear  waves 
with  nominal  resonance  of  2.25  MHz  ,  an  oscilloscope  (Tektronix 
model  455)  »  auxiliary  attenuators  ,  and  a  pneumatic  fixture 
that  provides  constant  pressure  between  transducers  and  specimen 
interfaces.  Couplants  AET  SC-6  and  Panametrics  SWC  were  used  at 
interface  of  transducers  and  specimens  for  longitudinal  and 
shear  tests  respectively. 

The  transmitting  transducer  is  excited  with  a  tone  burst 
voltage  signal  of  amplitude  100  Volts  (peak  to  peak);  this 
signal  imposes  a  stress  wave  (assumed  plane)  onto  the  specimen. 
After  travelling  through  the  specimen,  the  signal  is  captured  by 
the  receiving  transducer  and  displayed  on  the  oscilloscope. 

The  velocity  measurements  of  ultrasonic  waves  in  the 
through  transmission  configuration  may  be  affected  by  many 
parameters,  including  the  transducer-specimen  interface  pressure, 
the  mechanical  properties  and  thickness  of  the  couplant,  the 
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alignment  between  the  specimen  and  the  transducer, 
nonparallelism  of  the  faces  of  the  specimen  and  the  inherent  time 
delay  of  the  electrical  system.  To  avoid  these  effects,  the 
following  steps  were  taken: 

a)  The  specimen  was  very  carefully  centered  on  the  transducer 
face . 

b)  The  interface  pressure  was  maintained  at  4.2  Bar  during 
all  the  teats.  This  value  is  slightly  above  the  saturation 
pressure  [5]  and  guarantees  a  stable  output. 

c)  The  specimens  were  machined  within  established  tolerance 
levels  for  parallelism  between  faces  . 

d)  To  evaluate  the  velocity  ,the  time  shift  between  the  same 
individual  cycle  of  the  input  and  ouput  signals  was  measured.  A 
typical  example  of  the  input  and  ouput  signals  is  shown  in  Fig. 3. 
A  few  leading  cycles  were  not  used  because  of  their  transient 
nature.  In  order  to  cancel  the  inherent  time  delay  due  to  the 
electrical  system,  otherwise  identical  specimens  of  different 
lengths  were  used.  If  the  time  shifts  for  specimens  of  lengths 

are  t^,  respectively,  the  phase  velocity  v  can  be  expressed  as 

V  =  (L^  -  )/(t2  -  )  (1  ) 

The  tests  were  conducted  at  different  input  frequencies: 
0.28,  0.50  ,1.0,  1.5  and  2.0  MHz  .  The  selected  frequency 

range  depended  on  the  frequency  range  of  operation  of  the 
transducers  and  on  the  purpose  of  application  of  the  values 
measured . 
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Specimens 


The  material  used  was  3M  Scotchply  1002,  continuous  glass 
fiber  (type  E  )  in  a  165°C  (329°F)  curing  epoxy  matrix.  The 

composite  was  cured  in  a  heated  press  at  0.69  MN/m  (100  psi) 
according  3M  specifications  [12]  .  The  resulting  resin  content 
was  36  percent  by  weight.  Initially  a  25.4  mm  (1  in)  thick  plate 
with  dimensions  355.6mm  x  279.4  mm  (l4in  x12in)  was  fabricated. 
After  removing  the  edges  (approximately  10  mm  ),  the  specimens 
were  cut  from  the  plate  with  a  diamond  saw.  The  positions  of  the 
specimens  with  respect  to  the  longitudinal  axis  of  the  plate  (z) 
are  shown  in  Fig  .4  . 

The  cross-sectional  dimensions  were,  based  on  previous 
experience  [5,13],  smaller  than  the  diameter  of  the  transducers 
faces.  All  cross-sections  were  square  with  12.7  mm  (0.5  in) 
sides . 

The  length  to  be  traversed  by  the  waves,  ranged  from 
5mm  (0.19in)  to  29.6  mm  (l.06ln).  These  dimensions  were  limited 
by  the  geometry  of  the  original  plate  and  the  positions  of  the 
specimens  with  respect  to  the  longitudinal  axis  of  the  plate 
Attention  was  given  to  the  degree  of  parallelism  of  the  final 
surfaces  rather  than  to  the  dimensions  ,  since  the  last  could  be 
measured  a  posteriori  for  the  velocity  calculations.  The  degree 
of  parallelism  (along  the  length)  was  maintained  better  than 
+.013  mm  (+.0005in)  on  all  specimens. 

Since  the  velocities  are  to  be  measured  in  5  different 
directions  with  respect  to  the  longitudinal  axis  of  the  plate 
(0°,  30°,  45°»  60°,  90°),  at  least  10  specimens  (2  specimens 
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per  direction)  would  be  required  to  make  all  the  measurements 
However,  to  allow  for  possible  errors  and  variations  in  the 
material,  three  specimens  of  different  lengths  were  used  for  each 
of  the  directions. 
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RESULTS 


The  times  measured  for  the  waves  to  travel  through  a  certain 
thickness  of  material  were  found  to  be  independent  of  the 
frequency  for  each  of  the  investigated  directions.  This  finding 
shows  that  the  material  is  nondispersive  in  the  frequency  range 
used.  The  wave  speeds  were  calculated  from  the  time  measurements 
in  each  direction,  and  the  results  were  averaged.  The  confidence 
level  of  time  measurements  is  95  percent.  The  results  are  shown 
in  Table 

The  largest  values  for  the  velocities  occurred  for  the 
imposed  longitudinal  waves,  P  ,  followed  by  the  shear  waves  with 
polarization  in  the  xz  plane,  SV,  and  finally  by  the  shear  waves 
with  polarization  in  the  y  direction,  SH. 

For  the  purpose  of  comparison,  the  values  of  velocities 
along  the  principal  directions  were  also  calculated  from  layered 
media  theory  [14]  and  are  also  shown  in  Table  3*  Table  3  shows 
also  the  percentage  differences  between  the  experimental  and 
layered  media  theory  values. 

From  the  experimental  values  of  the  velocities  ,  the  elastic 

constants  of  the  stiffness  matrix  were  calculated  using  the 

expressions  in  Table  1 .  The  value  used  for  the  density  was 
3 

1850  kg/m  according  to  manufacturer's  data  [12].  The  final 
results  for  the  elastic  constants  of  the  stiffness  matrix  are 
shown  in  Table  4  with  the  corresponding  elastic  moduli  and 
Poisson's  ratios  shown  in  Table  5.  Observe  that  the 
measurements  along  nonprincipal  directions  of  propagation  (such 
as  the  directions  30°,  45°  and  60°  with  respect  to  z  )  are 
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supposed  to  give  values  of  velocities  that  can  be  used  for  the 
calculation  of  the  same  constant  C._  (see  Table  1).  The 

apparently  excessive  number  of  measurements  was  done  because  for 
direotions  other  than  the  principal  ones  (x,  y  and  z) ,  the  error 
in  the  time  measurements  can  be  large.  This  is  due  to  the  fact 
that  the  directions  of  plane  waves  originally  imposed  onto  the 
material  are  modified  by  a  deviation  effect  [15]  »  possibly 
introducing  errors  in  measured  times.  The  final  value  of  the 
constant  was  obtained  by  averaging  the  individual  values 

of  C._  calculated  for  each  of  the  nonprincipal  directions 

15 

(50  ,  45  .t  60  j  measured. 

In  general,  whenever  the  same  elastic  constant  of  the 
stiffness  matrix  could  be  calculated  from  measurements  along 
different  directions,  the  final  value  was  obtained  by  averaging 
the  values  obtained  for  that  constant  in  each  individual 
direction . 
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CONCLUSIONS 


The  measured  phase  velocities  indicate  that  the  fiberglass 
epoxy  material  tested  is  nondispersive  for  frequencies  between 
0.28  and  2.0  MHz. 

The  comparison  of  the  experimental  results  with  the  results 
of  layered  media  theory  show  reasonable  agreement  for  the  values 
of  velocities  of  longitudinal  waves.  For  shear  waves,  the 
percentage  differences  between  some  of  the  experimental  and 
theoretical  values  is  large. 

The  values  of  the  velocities  were  used  to  calculate  the 

elastic  constants  of  the  stiffness  matrix,  elastic  moduli  and  the 

Poisson’s  ratios  for  the  fiberglass  epoxy.  The  values  found  are 

elastically  compatible.  For  the  compatibility  of  elastic 

properties  the  Poisson’s  ratios  v  and  v  must  obey  the 

^  xy  zx 

following  condition  [16] 

-1  ^  V  6  1 

(2) 

V  i  (E  /E 
zx  '  x'  z' 

where  E  and  E  are  the  extensional  moduli  of  the  composite  in 
X  z 

the  X  and  z  directions,  respectively. 
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tepiesents  the  density  end  v  the  phase' velocity  in  the  propagation  direction 


TABLE  2  Relations  between  Elastic  Constants,  Elastic 
Moduli  and  Poisson’s  Ratios  for  Transversely 
Isotropic  Medium. 
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(*)  Characterizes  the  contraction  in  the  x  direction  when  tension 
is  applied  in  the  y  direction. 
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TABLE  3  Velocities  for  Unidirectional  Fiberglass  Epoxy 


1 

Direction  j  Imposed 

of  j  Wave 

Propagation!  Type 

1 

1 

1 

1  Velocity  (m/s) 

1 

1  Percentage 
j  Difference 

1  Experimental 

1 

1  Experimental j 

1  1 

1  1 

Layered 

Media  Theory 

z 

1  P 

1 

i  4692.8  j 

4344.7 

j  7.4 

30°  (♦) 

1  P 

i  4181.0  i 

mm 

i  - 

45° 

1  P 

1 

j  3558.7  j 

— 

i 

o 

o 

VO 

i  P 

1 

i  2940.9  i 

— 

I 

X 

i  p 

i  2391.9  i 

2300.8 

1  3.8 

z 

1  SH 

1 

i  1550.4  i 

1380.6 

1  10.96 

30° 

i  SH 

i  1521.1  i 

— 

1  - 

45° 

i  SH 

1 

j  1480.6  i 

1  1 

— 

1 

60° 

I  SH 

1 

j  1399.5  j 

— 

j 

X 

I  SH 

1 

i  1343.0  j 

1380.6 

i  -2.7 

z 

j  sv 

1 

i  1586.3  i 

1  1 

2835.0 

i  -78.72 

o 

o 

KV 

1  sv 

1 

i  1704.7  i 

— 

1  - 

45° 

j  sv 

1 

j  1836.1  j 

1  1 

— 

j  - 

60° 

i  sv 

1 

i  1904.2  j 

1  1 

- — 

1  - 

X 

i  sv 

1 

1  1506.3  ! 

1  1 

1380.6 

i  8.34 

(*)  All  angles  are  with  respect  to  z  axis 


TABLE  4  Values  of  Elastic  Constants  of  the  Stiffness  Matrix 


Elastic  Constant 

1 

1  Experimental  Value  (GN/m  ) 

1 

i  10.584 

I 

I 

I 

Sa 

i  4.098 

I 

I 

S3 

i  4.679 

I 

I 

I 

=33 

i  40.741 

I 

I 

I 

^44 

i  4.422 

I 

I 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

VO  1 

VO  1 

0  1 

1 

1 

1 

1 

i  3.243 

I 

I 

I 
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TABLE  5  Values  of  Elastic  Moduli  and  Poisson's  Ratios 


Elastic  Moduli 
and  Poisson  ratios 


Experimental  Values 


Extensional 

n 

GN/m^ 

Eztensional 

j 

1 

1 

8.785 

GN/m^ 

Poisson 

i 

0.354 

xy  1 

1 

1 

Poisson 

^  1 

XZ  1 

1 

0.319 

Poisson 

V  1 

ZX  1 

1 

0.074 

Shear 

! 

xy  1 

3.243 

GN/m^ 

Shear 

G  1 

XZ  1 

4.422 

GN/m^ 
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pulsed  Oscillator 
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Fig. 2  Schematic  of  experimental  system  for  through  transmission 


loo  volts 


iimmimiiui 

ramn 

■■■■ 


■ig.3  '  Typical  input  (upper)  and  output  (lower)  signals  in 
through  transmission  measurements  showing  reading  of 
time  t.  using  the  third  positive  peak  of  each 
signal . ^ 
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